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Abstract 

A perturbation theory for Massive Thirring Model (MTM) in radial quan- 
tization approach is developed. Investigation of the twisted sector in this the- 
ory allows us to calculate the vacuum expectation values of exponential fields 
(exp iaif (0)) of the sine-Gordon theory in first order over Massive Thirring Models 
coupling constant. It appears that the apparent difficulty in radial quantization 
of massive theories, namely the explicite "time" dependence of the Hamiltonian, 
may be successfully overcome. The result we have obtained agrees with the ex- 
act formula conjectured by Lukyanov and Zamolodchikov and coincides with the 
analogous calculations recently carried out in dual angular quantization approach 
by one of the authors. 



1 Introduction 

It is almost two decades the two dimensional exactly solvable models of QFT at- 
tract much attention because of their wide applications in the condenced matter 
physics and string theory. From the other hand, the knowlidge of the exact solu- 
tions of such nontrivial interacting theories provide us with better understanding 
of the general concepts of QFT. While the complete on-shell solution (i.e. the 
mass spectrum and exact s-matrix) for many of 2d Integrable QFT'es (IQFT) are 
well known, the construction of form-factors and correlation functions are under 
current intensive ivestigations now. 

The main subject of investigation in this paper is the most studied example 
of 2d IQFT the sine-Gordon model, which is defined by the action 

SsG = Jd'xS^^d,ipd'ip + 2ficosPip^ (1.1) 

where cp is a. real Bose field. The spectrum of this model includes the soliton, 
anti soliton and some number (depended on the coupling constant (3) of their 
bound states named breathers. The scattering in this model is factorized: the 
many-particle scattering process is reduced to the two-particle ones and the sets 
of two-momenta of incoming and outgoing particles are identical. It is well known 
since 1975 that the SG model is equivalent to the Massive Thirring model (MTM) 
with the action 



d^x Uij^dui^ -Mijtp-^ (^7>) (ijj.ij) I , (1.2) 
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where ?/), t/; are two component Dirac spinors. This equivalence assumes the iden- 
tification of the fundamental (anti) Fermions of the action ( |1.2| ) with the (anti) 
solitons of sine-Gordon model and the following relations among parameters and 
currents, estabhshed by Coleman 0: 

f = ^-l;r.*r., = -£.-a„^ (1.3) 

More recently Al.Zamolodchikov has obtained an exact relation between the 
soliton mass M and the parameter \i in the action ( |1.1|) 
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In this paper we consider the Vacuum Expectation Values (VEV) of exponen- 
tial fields in the sine-Gordon model 

Ga = {expiaip{0)) , (1.6) 

where the exponential fields are normalized by the condition 

^^ia^ix)^-ia^iy)\^^^ -. |x - yf'"' ttS \x - y\ ^ 0, (1.7) 

which emphasizes that the UV limit of this theory is governed by the c = 1 
conformal free bosonic field. 

For two special values of sine-Gordon coupling constant, namely for f3 ^ 
(semi classical limit) and /?^ = 1/2 (free fermion case), this function admits a 
direct calculation. The authors of have used these special cases to guess 
the following expression for the expectation values (|1.6|) for generic (3^ < 1 and 
\Re{a)\<l/{2f5) 
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In order to support the formula ( |1.8|) , some extra arguments, based on the 
reflection relations with Liouville reflection amplitude have been presented in the 
subsequent papers 0, P, but there is no rigorous mathematical proof yet. The 
article |§] provides another evidence supporting the Lukyanov-Zamolodchikov for- 
mula (|1.8| ), where its correctness has been checked in flrst order of MTM coupling 
constant g, by us of perturbation theory in angular quantization approach. 

In present paper we apply radial quantization to the same problem. The 
Hamiltonian of massive theories in the radial quantization approach has explicit 
time dependence [0. It appeared that this apparent difliculty can be overcome. 
We believe that such perturbative calculations substantially increase the confl- 
dence in reflection relations method as whole, which appears to be very powerful 
tool for investigation of 2d Conformal Field Theory (CFT) and IQFT 0, [|1^ . 



The paper is organized as follows. In section 2 we introduce the radial quan- 
tization of MTM. In section 3 we calculate the VEV ( |1.6|) at free fermion point 
g = 0. The calculation of VEV in the flrst order of perturbation theory is pre- 
sented in section 4. Here a special attention has been paid to regularization 
procedure of the product of local flelds at the coinciding points, which has some 
new features in comparison with the ordinary quantization in Cartesian coor- 
dinates. It appears that the Hankel-transform is a useful tool to carry out the 
calculations of section 4. The relevant mathematical details are presented in 
Appendix. 



2 Radial Quantization of the Massive Thirring 
Model 

In two dimensional space it is convenient to use the following representation of 
Dirac matrices 



-t\ , . 701 



and denote the components of Dirac spinors as 



(2.1) 
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In this notations the action (|1.7|) in Euclidean space takes the factorized form 
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^lidipR + ipidipL - -J- (V'tV'R - i^ni^- 



+^Vt^L4^J , (2.3) 



where z = x'^ + ix^, z = x^ — ix^ are the complex coordinates on the Euclidean 
plane, d = d/dz, d = d/dz and (Pz = 2dx^dx'^ is the volume element. 

As we are interested in the VEV's of local fields /e'"'^^''-* V which have rotational 
symmetry, it is natural to use the polar coordinates 77, 6 defined by 



^'n+iB 
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(2.4) 

and interpret r^, 9 as Euclidean time and space coordinates respectively. 

Since the conformal dimensions of the Fermi fields ipL and ipR are (1/2, 0) and 
(0, 1/2), they behave under the conformal transformations ( |2.4| ) as 
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where C,=Ti + i6,^ = T] — i6. The same transformation lows hold for the fields 



Thus in {rj, 6) coordinates the action ( p. 3D becomes 
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Treating the radial coordinate 77 as a " time" , we get the Hamiltonian 

+2^Vl V-L^li^i?] • (2.7) 
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H = fde 





The usual canonical quantization scheme will bring us to the standard "equal 
time" anti-commutation relations 



{^|JL (e) , ^I {6')} = 6{e-e'), {^^ (e) , v^J^ (o')} = Ho-o')- (2.8) 

As usual, in order to develop perturbation theory one first has to solve the 
problem with unperturbed Hamiltonian (i.e. to ignore the last quartic term in 
(|2.7|) ). We found it easier to handle this problem in Schrodinger picture, instead 
of more conventional in QFT Heisenberg or Interaction pictures. Thus our field 
operators ipL,R will not depend on "time" r] and the state vectors will evolve 
according to the Schrodiger equation. Let us define the creation, annihilation 
operators c\, dl, Ck, dk through the Fourier mode decompositions 
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where all sums are taken over all positive half-integers [M is the set of positive 
integers) . 

As a consequence of equations ( |2.8| ) and ( p.9| ) , one can easily get the following 
anti-commutation relations for the operators c^, (ifc, c|., (i|.: 

{cfc, ci} = {4, 4} = {4, rfj = {4> 4} = 0, 

{cfc,4}=4,«, {4,4} = '^fc,h k,leZ--, (2.10) 

where Z is the set of integers. As usual the Fock space (let us denote it H) 
has the following basic vectors 



n (4)" (4)" 10), (2.11) 

where n^ G {0,1} {hk G {0,1}) are the occupation numbers of "particles" 
created by the operators c^ {d\) out of bare vacuum |0), which by definition 
satisfies the conditions 

Cfc|0) = 4|0) = 0; keZ-]^. (2.12) 

In terms of the creation, annihilation operators is the Hamiltonian (|2.7|) ac- 
quires the form H = Hq + Hint, where the quadratic part Hq is equal to 

^0 = Yl [k (4cfc - 44 + clfcC-fc - d-kdlk) (2-13) 

-iMe"^ (cI4fe - d^kCk + dkC-k - c^-kd\ 
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The evolution of arbitrary state \s) along Euclidean time r] caused by the 
Hamiltonian Hq is given by the Schrodinger equation 

d 

-r—\s,r) = Ho\s,r). (2.14) 

Here and henceforth we prefer to use r = Me'' rather than rj. To find the 
general solution to the Schrodinger equation (|2.14 ) let us denote, that the Hamil- 
tonian Hq has a factorized form 

Ho= E K^ + i/f), (2.15) 

where the operator Hi {Hi ') includes only Ck, d-k, cl, d_i^ (4, c_fc, 4, c_^). 
This makes convenient to represent the Fock space Ti as an infinite tensor product 

,(1) ^^(2) 



n = 0kex-i/2[K'(»n):'), (2.16) 

where T-lj, and Hi. are four dimensional vector spaces with base vectors 

|0i'^ 44j0i^)), (even sector) 

4|0i'^), 4,|0«), (odd sector) (2.17) 



and 



lof), 44j0l')), (even sector) 

4|0f^), 4,|0l')), (odd sector) (2.18) 



respectively. The vectors |0^ ) and |0^ ) are defined by the conditions 
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for any k G M — 1/2. Note that the bare vacuum introduced earher (see (|2.12| ) 
) is equal to 



|0) = ®fc6Ar- 
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The operator iJ^ (//^ ) nontrivially acts only on the factor 1-L\. {Hi ) of 
the full Fock space H (|2.16|) , hence we have reduced the initial QFT problem of 
infinitely many degrees of freedom to the simple quantum mechanical one, with 
four dimensional Hilbert space. A further simplification provides the observation, 
that the reduced Hamiltonians H^ , H^ don't mix even and odd sectors (see 



(|2.17| ), ( p.lSp ). The resulting Schrodinger equations in this reduced spaces take 
the form 
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a,ir) + P,ir)cidU] |0«) = //f^ [a,(r) + /3,(r)cldlj |0«) = 

'-kak{r)+trPkir) + {k(3u{r) - irak{r))c{d}_j\ |oi'^), (2.21) 



lk{r)c\ + 5,{r)dl,] |0«) = i7« k(r)cl + 5,{r)dl,] lO^) = 0. (2.22) 



Evidently, to obtain the equations for the another sector with Hamiltonian 
Hi. , one simply has to change the upper indices (1) into (2) and make the 
substitutions c^ — > d^ and d_k -^ c_fc. Thus, the differential equations for the 
unknown functions afc(r), Pk{r), 'yk{r), 5k{r) in both cases remain the same: 
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The second pare of these equations show, that in fact 7jt and 5^ are constants, 
while the first pare reduces to the modified Bessel differential equation with 
general solution 
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One should fix the constants 0^,6^,7^ and 6k imposing initial conditions at 
the arbitrary "time" tq. For the further application let us write down explicit 
expressions with specified constants for two basic cases 
a) when the initial state coincides with |0^ ) or |0^ ) : 

ftfc (r) = ^/F¥^ (i^fc+i (ro) 4_i (r) + 4_^i (ro) Ki^_i (r)) , 

Pk (r) = tV^o {Kk+i (ro) h+i (r) - 4+| (ro) i^^+i (r)) , (2.25) 



(t)^(t)in(i)\ 
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b) when the initial state coincides with ^c^ w_^|0^ ) or ^rf^ c_^|0^ ) 

ak (r) = y/¥¥^ (k^^_^_ (ro) 4_i (r) - 4_i (ro) K^^^ (r)) , 

/?fc (r) = i^¥F, {k,_._ (ro) 4+1 (r) + 4_ 1 (ro) iT.^i (r)) . (2.26) 



During the proof of the formulae ( p.25|) ( p^26D we have the standard Wronskian 
identity for the modified Bessel functions 



4-1 (r) i^,+i (r) + 4+1 (r) i^fc 1 (r) 
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(2.27) 



It is interesting to note, that due to explicit time dependence of the Hamil- 
tonian, the system being initially at the ground state of that particular moment, 
after finite time of evolution will find himself in an excited state. Nevertheless, 
long time evolution of any state with non vanishing overlap with the ground state 
of the initial time, eventually approaches to the ground state of the infinite future 
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Evidently, at the small r's (far past), the ground state approaches to the bare 
vacuum |0). In particular, if r ^ 1 and ro -C 1 (|2.25|) gives 
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3 The VEV's of the Exponential Fields in 
Free Fermion Case 



The VEV (|T_, 



with SsG being the action (|1.1| ), can be expressed alternatively in terms of the 
appropriately regularized (see below) Euclidean functional integral over the Dirac 
fermions 



p—-A.MTM 



Gia) = ^ -. (3.31) 



where Amtm is the Euclidean action (|2.6| ). The functional integral in the nu- 
merator of (|3.31| ) is taken over the space Ta of the twisted field configurations 
■0 (^, ^) and ^ (2;,^), which acquire the phase 
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^ (z, Y) -^ e~^ {z, z) , ?/> (z, z) -^ e~'~ilj (z, z) , (3.32) 

when continued analytically around the point z = in anti-clockwise direction 
[0. This is due to the non-trivial monodromy of Dirac fields with respect to 
the exponential fields expiaip (0) . It is easy to see, that to ensure such twisted 
boundary conditions on Dirac fields, one has to shift Fourier mode indices as 
follows 

k ^ k — —, in sector 1 (i.e. in Ck, dk sector) 

P 

k —^ k + —, in sector 2 (i.e. in dk, Ck sector). (3.33) 

P 

For example, the Fourier decomposition of the field ^l(^) takes the form (cf. 
with the first line of (g^)) 



,P,{9) = ^ Y: (c._„ e-^^'^-")^ + 4_„ e^^'^-")^) , (3.34) 

V2vr ^^^_, 

where 

a^^. (3.35) 

With such shifts, all the results of the previous sector remain valid since we have 
never used the arithmetical properties of the Fourier mode indices. 

In radial Hamiltonian formalism the regularized version of the functional in- 
tegral (|3.31|) may be represented as 

{oo\S{r,ro) |0) 
8 



where the matrix element of the evolution operator S{r, tq) in the numerator is 
taken in twisted sector (this is indicated by the lower index a) . To regularize the 



expression, in ( 3.36 ) we have assumed, that the evolution begins at some small 
ro- A simple Conformal Field Theory consideration^, which takes into account 
the fact, that the conformal dimension of the field e*""^ is a^, leads to 



G.(a) 



lim (ro) 

ro— >0 



-2a^ 



Go(a,ro) 



(3.37) 



In general case it is not known yet how to calculate the functional integral ( |3.31| ) 
or the matrix elements in (|3.3(j| ) exactly. Below we'll evaluate ( ^.36 ) at the free 
fermion point g = 0. As in this case we already know the time evolution of every 
constituent of the vacuum |0) (see eq. ( p.20| )) from the previous section, it is not 



difficult to pick up all the necessary factors from (|2.29|) with appropriate shifts 
of Fourier mode indices and obtain 



Go{a,ro)= (3.38) 
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where we have added a subscript to G in order to emphasize, that the free 
Fermion case g = is considered. We have to be careful, when evaluating infinite 



products in ( p.39| ) and treat the ill defined sums like J2{'i ^ o-) by means of 

i=0 

Riemann (^-function regularization. Let us remind that 



C{z,a) = Y, 



i=0 



i + a) 



and 



(3.39) 



C (-1, a) + C (-1, -a) - 2C (-1, 0) = -a^ (3.40) 

To carry out the remaining infinite products of F-functions (also divergent, if 
treated literally), it is convenient to use the integral representation 
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As a result we obtain simple geometric progressions, coming from the first term 
of the eq. ( |3.41|) and contributions, coming from the second term, which can be 



^In the limit r — > the action (2^) describes the Massless Thiring Model, which is well 
known to be conformal invariant. 



easily handled applying C- function regularization once more. The final expression 
has the form 



Go(a,ro) 



^0 



exp 



sinh^(at) _ ^ _2t 
sinh t 



dt 
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(3.42) 



or, taking into account equations (|3.37|) and ( p.35|) with the free Fermion point 
value p = I/V2 



Gaia) = i y) 
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(3.43) 



This is in full agreement with the result, obtained by S.Lukyanov and A.Zamolodchikov 
in 0], using angular quantization technic. 



4 VEV of the Exponential Field in the First Or- 
der of Perturbation Theory 

In this section we calculate the VEV ( |1.6| ) in first order of the MTM's coupling 
constant g. The perturbation is given by the last term of the Hamiltonian 

27r 



n,nt = 2gjN{^t^L'^j,^R)de, 



(4.1) 



where we denoted by N{. . .) an appropriately regularized product of local fields 
at coinciding point. One has to chose such a regularization, which will not break 
the translational invariance of the theory if it is transformed back to the initial 



Euclidean coordinates x^, 



X 



The conventional normal ordering with respect 



to creation-annihilation operators fails to satisfy this condition, because of the 
non-trivial time dependence of physical vacuum. Instead, we should suppress all 
the contractions among fields inside the correct normal ordering symbol N[...) 



■.^i^L\{:^UR\, 



+ 
(4.2) 



where the vacuum expectation value of any operator X is defined by 

(oo|5(i?,r)XS(r,ro|0)) 



(^)c 
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(4.3) 
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with all matrix elements taken in untwisted sector. The first part of eq. 
could be understood for example as a zero distance limit of corresponding point- 
split expression. In ( |4.3|) a small initial time tq and a large final time Rq are 
introduced in order to keep intermediate expressions finite. R and r^ eventually 
will be sent to and cxd correspondingly. Note also the explicit r dependence of 

which reflects the inhomogeneity of "time" in our scheme 



([4.3| ) and hence of (^ 

of quantization. The standard time dependent perturbation theory in first order 

of the coupling constant g gives 



G'(o,ro)= lim | (oo|5 (i?,ro) |0)^ + / (oo|5 (/?,r) ff,„,5 (r,ro) |0), 



ro 



(4.4) 

As we already have obtained explicit expressions for time evolution of states from 
various sectors of Fock space in section 2, it is not difficult to calculate the matrix 
element under the integral in eq. ( [4.4| ) 



G(a,ro)= lim (oo |5 (i?, ro)| 0), {1- 



R 



^ / rdr [2Ik+l-aKk-ah+l+aKl+a — h+l-aKk-ah+l-aKl-a — 



Ik+1+aKk+ah+l+aKlJ^a — h-aKk-ah-a+lKl^a+l — h+aKk+ah+a+lKl^a+l — 
Ik+l-aKk+l-ah+l+aKl+l+a — Ik+aKk+ah-aKl-a + -^fe+l-^fe+l-^Z-a+l-^i-a+l + 
hKkk-a+lKl-a+l + hKkk+aKl+a + hKkk+a+lKl+a+l + hKkk-aKl^a + 
Ik+lKk+lIl^a+lKl^a+l + h+lKk+lh-aKl-a + h+lKk+lh+aKl+a — 

h+iKk+ik+iKi+i - Ik+iKk+ikKi - IkKkk+iKi+i - IkKkkKi]} , (4.5) 



with the pre factor (oo|S'(/2, ro) |0)a given by eq. (|3.42| ) (in (|3.42| ) we have to insert 
(3 = 45(1 ~ ^ + o{g)) and expand the resulting expression over g up to linear 
term). The remaining calculation of integrals over the quartic Bessel functions is 
presented in the appendix. Using these results we obtain 
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(4.6) 



X 
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Now performing summation over k and / in the third line of the eq. ([4.6|) we obtain 
logarithmically diverging at t = integral. In fact, the same problem we have 
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faced when carrying out calculations exactly at the free-fermion point. Indeed the 
the product in ( p.39| ) diverges for large k, I, but we overcame this difficulty using 
C-function regularization inside the integral representation of F-function (3.41). 



Here we'll not care of a similar appropriate regularization. Rather, noticing that 
various regularization scheme will differ from each other by a term ~ a^, and that 
the coefficient of — l/2a^ in the expansion of (e*'*'^) is just the VEV (0^), which 
has been calculated in 0] using standard Feynman diagram technic with result 
(below 7 = 0.577216... is the Euler constant) 



,2^r.^^ ./.,., 1_„.,./r.NN , 9 



(0)) = -4(1 + 7 + log(M/2)) + ^(7C(3) - 2) + Oig'), (4.7) 



n 



we simply cut the above mentioned integral over t on the lower bound and equate 
the undefined coefficient of — l/2a^ to the one predicted by the eq.([4.7|). The final 
result is 




' exp ^ / ( ^ X - " e"'* ]^} y< (4.8) 



M\°' \ ffsmh\at) ^2„-2t\ dt 
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asinh(2at) sinh^(at) 



2 / W \ sinh t It 







2 sinh^ t sinh^ t 



dt- 2a' \og2 



+ 0{g' 



which is in complete agreement with the Lukyanov-Zamolodchikov formula ( |1.8| ). 



Acknowledgments 

We are grateful to A.Sedrakyan for useful discussions. R.P. thanks theory division 
of LAPP and especially P.Sorba for warm hospitality and discussions. The work 
of V.M. and R.P. was partially supported by INTAS grant 96-690. 



Appendix 

It appears that the Hankel-transforms are appropriate tools allowing us to per- 
form the integration over r in (|4.5| ). Roughly speaking, in polar coordinates the 
Hankel-transforms play the same role, as the ordinary Fourier-transforms in the 
Cartesian one. 

Let us briefly recall the main formulae concerning to the Hankel-transforms 
(for details see [|13| and references therein). The z/-th order [u > —1 ) direct and 
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inverse Hankel-transforms of the function f{x) defined on (0, oo) are given by 



f (x) = Ju (sx) fu (s) sds, (A.l) 



oo 

fu (s) = I Ju {sx) f (x) xdx, (A. 2) 



where Jj, is the Bessel function. In complete analogy with the case of Fourier- 
transform, it follows from ( [A .11 ), ( [A .21 ), that the "scalar product" of any two 



functions f{x), g{x) coincides with that of their images 

3 OO 

/ {x) g {x) xdx = fu (s) gl (s) sds. (A. 3) 



To use (|A.3|) for the calculation of the integral in eq. ([4.5|) we need to know 
Hankel images of the functions Iy{x)Ku{x) and Iy^i{x)Ky[x) which can be easily 
obtained from the general formula [|13| 



K^y ix) la ix) 



-V 



3 



/'j_^+^(2xsinht)e-('^+^)*rft, 



-i?e(z/ + /i) < -, i?e(-z/ + /i) > -1, (A.4) 



namely 



oo 

h {x) Ki (x) = fjo (xs) ^^=e-2'*(«)srfs, (A.5) 

oo 

Ii+, (x) K, (x) = fj, (xs) _==e-(2'+i)*(^)srfs, (A.6) 

where t{s) is defined by 

ds 
2 sinht = s, dt = , (A.7) 

Though the direct application of eq. ( |A.3| ) to each term of ( [4.5[ ) at first sight 
seems to be problematic due to the logarithmic divergence of the integral over at 
large r, but nevertheless it leads to a correct result, because of mutual cancellation 
of these divergences by various terms. 

^Sincc the functions we are dealing with are regular in the interval (0,cxd), the only thing 
one has to care of is the convergence of integrals at the extreme points and cx). 
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